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ABSTRACT 

Simulations of generic pushbroom satellite 
hyperspectral sensors have been performed to 
evaluate the potential performance and validation 
techniques for satellite systems such as COIS 
(NEMO), Warfighter- 1 (Orb View-4), and Hyperion 
(EO-1). The simulations start with a generation of 
synthetic scenes from material maps of studied 
terrain. Scene- reflected radiance is corrected for 
atmospheric effects and convolved with sensor 
spectral response using MODTRAN 4 radiance and 
transmission calculations. Scene images are further 
convolved with point spread functions derived from 
Optical Transfer Functions (OTF s) of the sensor 
system. Photon noise and detector/electronics noise 
are added to the simulated images, which are also 
Finally quantized to the sensor bit resolution. Studied 
scenes include bridges and straight roads used for 
evaluation of sensor spatial resolution, as well as 
fields of minerals, vegetation, and manmade 
materials used for evaluation of sensor radiometric 
response and sensitivity. The scenes are simulated 
with various seasons and weather conditions. Signal - 
to-noise ratios and expected performance are 
estimated for typical satellite system specifications 
and arc discussed for all the scenes. 

Key words: hyperspectral remote sensing, sensor 
design, computer simulations, virtual prototyping 

/, INTRODUCTION 

Hyperspectral imaging is an exciting, emerging, 
remote sensing technology that enables one 
instrument to have utility for a wide range of 
applications. Because hyperspectral sensors record 
hundreds of spectral bands for each location they 
image, the data they provide significantly improve 
capabilities for automated atmospheric correction, 
spectral unmixing, and feature extraction, which form 


the basis for many remote sensing applications. In the 
near future, three satellites with hyperspectral sensors 
are going to be launched into low' earth orbits. The 
instruments include the Coastal Ocean Imaging 
Spectrometer (COIS) on the Naval EarthMap 
Observer (NEMO) satellite. Warfighter- 1 on the 
Orb View-4 commercial satellite, and Hyperion on the 
NASA's Earth Observing- 1 (EO-1) spacecraft. 
Despite some distinct differences between the 
instruments, their ability to provide hyperspectral 
image cubes from space-borne platforms is the 
common feature which encourages comparison and 
cross-validation between them. AVIRIS, an airborne 
hyperspectral sensor used on the ER-2 and Twin 
Otter platforms, is a demonstration of the outstanding 
radiometric performance, which can be achieved by 
such a system (Green cl al. 1998). To understand 
w r hal kind of performance one can expect from the 
satellite hyperspectral sensors, we have created a 
computational prototype of a generic hyperspectral 
sensor on a satellite platform. In recent years, rapid 
advances in computer technology have allowed 
computational modeling and simulations to become 
valuable tools in predicting performance 
characteristics of future instruments and devices 
before physical parts are made and assembled 
(Rlonski et al. 1997). The Virtual Product Laboratory 
(VPL) software environment developed at the NASA 
Stennis Space Center allows for such studies of 
remote sensing systems (Gasser et al. 1999), 

//, INSTRUMENT DESIGN 

Our virtual hyperspectral instrument is a pushbroom 
system with a reflective telescope, a prism 
spectrograph, and a state-of-the-art focal plane array 
detector optimized for the visible and SWIR spectral 
regions. The sensor is designed to provide imagery 
with a 30- m ground sampling distance (GSD) from 
the altitude of 705 km. Because the focal plane array 
is supposed to have detector pitch of 24 pm (with 
100% fill factor), the required effective focal length 
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of the optics is 564 mm. Therefore, the instantaneous 
field of view (IFOV), defined here as a sampling 
interval not as the full width at half maximum 
(FWHM) of a point- spread function, is 0*0426 mrad. 
Additionally, because we also assume that the 
detector array has 640 pixels in the spatial (cross- 
track) direction, field of view of the sensor is 1 .5°. A 
reflective telescope with the Cassegrain design and 
173 is the most appropriate for such a case (Smith 
1992). This tY# implies aperture size of 188 mm. 
Optical efficiency of the combined telescope - 
spectrograph system is conservatively assumed to be 
10% over the entire spectral range. The optics are 
supposed to be only diffraction- limited, with no 
aberration blur. Spectral range of the sensor is from 
400 nm to 2500 nm, with the sampling interval of 10 
mn -211 bands with the Gaussian shape and FWHM 
of 10 nm. The detector array and electronics are also 
characterized by an integration time of 4 ms and a 
12-bit quantization. A simple noise model, which 
accounts for the detector and electronics noise (the 
illumination- independent noise floor) as well as the 
photon (shot) noise, is applied in the simulations. The 
noise floor is assumed to be 5 nW / (cm 2 sr nm). The 
satellite is additionally supposed to introduce random 
jitter with an RMS blur of 0.0242 mrad* Although 
without much detailed knowledge of a specific 
design, and without simulations on a sub-system 
level, it is difficult lo estimate the actual performance 
parameters such as the optical efficiency or the noise 
floor, we have selected values, which we believe 
would be typical for present-day hyperspectral sensor 
systems. Table I lists the parameters of the simulated 
sensor* 


Table L Characteristics of the simulated satellite 
hyperspectral sensor system 


Altitude 

705 km 

GSD 

30 m 

Spectral range 

400 - 2500 nm 

Spectral sampling 

10 nm 

Spectral resolution 

10 nm 

Number of spectral bands 

211 

Spatial sampling (IFOV) 

0.0426 mrad 

Field of view' 

1,5° 

Detector pitch 

24 pm * 

Effective focal length 

564 mm 

m 

3 

Aperture diameter 

1 88 mm 

Optical efficiency 

10% 

Integration time 

4 ms 

Quantization 

1 2 hits 

Noise floor 

5 nW / (cm 2 sr nm) 

Random jitter RMS hlur 

0.0242 mrad 


IlL SENSOR MODEL 

A system level approach based on Optical Transfer 
Functions (OTF's) is applied in the modeling of the 
sensor (Kopeika 1998). System OTF is formed by 
multiplication of individual functions taking into 
account the following factors: 

* Detector size 

* Diffraction 

* Linear motion 

* Random jitter 

* Atmospheric turbulence 

The system OTF is assumed to be separable in two 
orthogonal axes coincident with the detector array 
axes, one of whieh is parallel {and the other is 
perpendicular) to the direction of the satellite ground 
track. In both cross-track and along- track directions, 
the individual OTF’s are described by formulae 
presented below. Using detectors of finite size to 
spatially sample a scene results in the detector OTF is 
given by the equation (Holst 1998): 

OTFWutM = sin z(md 1 f ) 

where u is the object- space spatial frequency, d is the 
detector size, and /is the effective focal length of the 
optics. The sine function is defined here by: 
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Diffraction in the optical system results in image 
blurring described by the following function (Holst 


1998): 
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with Wo = D f A, where D is the aperture size and A is 
the wavelength. OTF associated with the orbital 
motion of the satellite is given by (Holst 1998, 
Kopeika 1998): 

= sinc(mn’T/ h) 


where h is the altitude of the satellite above the 
earth’s surface, v is the speed of the subsatellite point 
on the ground* and r is the integration time. The 
motion OTF affects the system OTF only in the 


along-track direction. The random jitter OTF is 
approximated with a Gaussian function (Holst 1998): 

OTF; tt >)=exp(-2/rW) 

where b is an object-space RMS diameter of the blur 
created by the random motion. The expression for the 
atmospheric turbulence OTF is based on the Fried’s 
approach (Fried 1966, Bel and 1993): 

0 T F,U-t,ul C n,c<' r ') 

= exp{- [f n f )f'Ou / r 0 ) i,3 [l - x-(Au / D) l/3 ]} 
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where r„ is called the atmospheric coherence length 
or the Fried’s seeing parameter, F is the Gamma 
function, and the factor x'is specified in the following 
way: 

• K— 0 for long exposure times (T> 10 ms) 

• K — () r 5 for short exposure times (r< 10 ms) 
and far- field propagation (h > D 2 !A) 

• K= I for short exposure times (r< 10 ms) 
and near-field propagation ( h < D'fX) 

Note: [MF(f )f /6 ~3,44. In this study, a value of 20 

cm, which is better than a typical one, is assumed for 
the Fried’s parameter (Behind 1993). 

Figures 1 and 2 show absolute values of OTF’s 
calculated for the simulated sensor. 
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Figure 1. Cross-track modulation transfer 
function of the simulated sensor for hand 16 
(550 nm). 


Figure 2. Along-t rack modulation transfer 
function of the simulated sensor for hand !6 
(550 nm). 


IV : IMAGING SIMULATION 

Both cross-track and along-track components of the 
system OTF arc used Lo derive a point spread 
function (PSF) for the sensor using an inverse Fourier 
transform. The calculations are performed in a 
discrete fashion using the fast Fourier transform 
(FFT) algorithm and a series of spatial frequencies 
dependent on ground sampling distance with which 
an input scene is defined. Range of the PSF is 
truncated to a rectangular area where the cross- track 
and along-track components of PSF are no less ihen 
1% of their maximum values. To protect radiometric 
accuracy of the simulations, the truncated PSF is 
normalized in such a way that a sum of all the 
elements of the PSF matrix is equal to one. The input 
scene image is convolved with the PSF to generate an 
image created by the simulated sensor. The 
convolution is performed in the spatial domain 
because wc found that this approach offers better 
performance than the one based on two-dimensional 
FFT. The savings are created when the convolution is 
performed as a block operation with the PSF kernel 
moving over the scene image with step defined by the 
sensor sampling distance. To ensure quality of the 
simulations, the input scene must be created with 
spatial sampling significantly finer than that of the 
sensor. Recent studies have shown that simulations 
with the sampling ratios of 3 to 7 produce adequate 
results (Jacobs & Edwards 1999). However, only 
sampling in the direction perpendicular to an edge 
was considered in that work. Our tests have indicated 
that much larger oversampling is needed to reliably 
simulate images of edges oblique to the sampling 
direction. Therefore, an oversampling of 15 was used 
in the present simulations. However, one must realize 


that such an oversampling comes with a significant 
computational price, especially in the amount of disk 
space required to store the scene image cubes. There 
is definitely a need for further studies on how to 
achieve the desired fidelity of the simulations while 
keeping the computational requirements in check. 
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Figure 3. Paint spread function of the simulated 
sensor for band 16 ( 550 run). The mesh is 
drawn using the sampling interval of an input 
scene. 


square root of the number of photons. Finally, the 
number of photons is converted back into radiance by 
reversing the above formula. Next, adding a random 
number generated according to the normal 
distribution with standard deviation equal to the noise 
floor simulates the detector and electronics noise. 
Intensity of the simulated image is additionally 
quantized to the number of levels given by the 
number of bits provided by the sensor for each pixel. 
The OTF and PSF calculations, the convolution, and 
the noise simulations are performed separately for 
each spectral band of the sensor, but at the end the 
output images are combined into an image cube, 

V SYNTHETIC SCENE 

While the scene radiance is convolved with the 
sensor spatial response during the imaging 
simulation phase, its convolution with the sensor 
spectral response is done during creation of the 
synthetic scene. In general, the simulation process 
(without effects of noise) can be summarized in the 
following formula (Schowengerdt 1997); 

l(x 0 ,y o,/U 

— JJJ L(x, y, X)R(X — Xn)PSF(x — x 0 , y — y a , X)dxdydX 


After the image convolution, detector and electronics 
noise as well as photon noise is added to each pixel 
of the simulated image. To simulate the photon noise, 
first the pixel radiance L is converted into the number 
of photons n using the following formula: 


n = 
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where / is the simulated sensor- generated image, L is 
the synthetic scene spectral radiance, R is the sensor 
spectral response, and PSF is the sensor spatial 
response (the point spread function). Assuming that 
the PSF does not change significantly for the 
wavelengths w r ithin a single spectral band, the spatial 
convolution can he separated from the spectral 
convolution to result in: 


where: 

£ 

optical efficiency 

S 

spatial sampling interval of an input scene 

^mas 

maximum value of the PSF 

D 

sensor aperture diameter 

H 

satellite altitude 

M 

spectral band width 

4 

spectral band center wavelength 

T 

integration time 

h 

Planck’s constant 

c 

speed of light 


Then, a random number, which is generated 
according to the Poisson distribution with the mean 
value equal to the number of photons, replaces the 
original value. For large numbers of photons, the 
Poisson distribution is approximated by the normal 
distribution with the mean value equal to the number 
of photons and the standard deviation equal to the 


^ (*o * So ’ ) 

= Jj L(x, .4,, )PSF{x - x (l , y - v 0 , X [t )dxdy 

w'here the in-band radiance is given by: 

L(x, y, 4 ) = J L(x, y , X)R(X - X 0 )dX 

Although, in our implementation, it is done in a 
discrete fashion, nevertheless, the former integral is 
calculated during the imaging simulations, while Lhe 
latter is evaluated during atmospheric radiative 
transfer calculations of the synthetic scene creation 
process. These calculations also incorporate into the 
simulations the geometric properties of the satellite 
motion: altitude above the Earth’s surface and 
azimuth of the ground track as well as orientation of 
the line-of-sight 


Image cubes of the synthetic scenes used as input for 
the simulations can be generated in one of the 
following three ways. The first approach uses the 
DIRSIG software package developed at Rochester 
Institute of Technology (Schott et a!. 1999), DIRSIG 
uses facetized models of terrain and objects such as 
buildings or vehicles to calculate scene radiance by 
performing ray tracing and atmospheric radiative 
transfer modeling based on MODTRAN. The second 
option is to use GCI Toolkit, which is a commercial 
software package developed by Photon Research 
Associates. GCI Toolkit uses rasterized models of 
terrain and clouds as well as atmospheric modeling 
based on MOSART. Availability of the rasterized 
terrain databases with sufficiently high spatial 
resolution limits our ability to use GCI Toolkit in the 
current simulations. The third approach is an ad hoc 
combination of a two-dimensional material map and 
the atmospheric radiative transfer model of 
MODTRAN, Such calculations are described in 
detail in another paper in these proceedings (Cao et 
ah 1999). 

For the si gnal-io- noise ratio (S/N) studies presented 
in this paper, flat uniform scenes, with the same 
material defined for all the pixels, w'ere created using 
the third method, MODTRAN 4 was executed using 
the correlated-k option (Berk et ah 1998) for the 
geographic location of the scene in the proximity of 
New' Orleans, Louisiana, at 30° N latitude and 90° W 
longitude. Because of the location, the mid-latitude 
summer model atmosphere was assumed for all the 
seasons, and as a result the only difference between 
the seasons is in the Earth to sun distance and the 
sun’s location in the sky. These differences are 
introduced by performing the calculations for the day 
of the year number 80 (spring), 172 (summer), and 
356 (winter). Other parameters included: C0 2 mixing 
ratio of 360 ppm v , water column scaling factor of 0.5, 
the default vertical ozone column, rural visibility of 
23 km, and a nadir looking from the top of 
atmosphere. Reflectance spectra of surfaces used in 
the simulations were obtained from the ASTER 
Spectral Library (version 1,1), w r hieh includes data 
from three sources: the Johns Hopkins University 
(JHU) Spectral Library, the Jet Propulsion 
Laboratory (JPL) Spectral Library, and the United 
States Geological Survey (USGS - Restoif) Spectral 
Library (for reference, see the WWW page at URL 
h ttp : // spec lib.jpl.na sa . go v ) . 

VL RESULTS 

The synthetic scene of the Rochester area was used to 
test capabilities of the hyperspectral imaging 
simulations. 



Figure 4. Image of a synthetic scene created 
using DIRSIG with I m GSD and the 21 1 
spectral hands of the simulated sensor. Bands 
23 (620 nm), 16 ( 550 nm), and 7 (460 nm) are 
displayed as the RGB colors, respectively , 
Prof John Schott Rochester Institute of 
Technology, provided the model geometry for 
this scene based on the Genesee River Gorge 
area in Rochester, New York. 

An image cube of the input scene consisted of 1024 
by 1024 pixels and 21 1 bands. Simulations of the 
imaging process do not change the number of 
spectral bands, while reducing the spatial resolution 
(see Figure 5). Effects of spatial sampling as well as 
optics and atmosphere are clearly visible when the 
simulated image (Figure 6) is compared with an 
image of the input scene (Figure 4). 



Figure 5. Image cube generated by the 
simulated hyperspectral sensor from the 
synthetic scene of the Rochester area 
created with DIRSIG, 



Figure 6. Image generated by the 
simulated sensor from the synthetic 
scene of the Rochester area created 
with DIRSIG . Bands 23 (620 nm) t 16 
(550 rim | and 7 f 460 nm) are displayed 
as the RGB colors f respectively. 


Image cubes created in the simulations performed 
with the uniform scenes for the S/N studies consisted 
of 100 by 100 pixels and 21 1 bands. For each of the 
bands, mean and standard deviation of image 
intensity were calculated using all the 10,000 pixels. 
The mean divided by the standard deviation was used 
as an estimate of the signal -to-noise ratio for the 
material used in the given uniform scene. Surfaces of 
the following materials were used in the sLudies: 

* water 

* minerals; 

al unite, cal cite, dolomite, kaolinite, quartz 

* soils: 

black loam, brown loam, dark loam, dune 
sand, gray clay 

* vegetation: 

dry long grass, lawn grass, sage brush, fir 
tree, pine tree 

* man-made materials; 

aluminum, asphalt, concrete, steel, terra 
cotta 

For each material, three scenes w r ere created w r ith 
atmospheric effects calculated using solar conditions 
for different seasons: w inter, spring, and summer. 
Spectral radiance and noise equivalent spectral 
radiance (MESR) obtained in the simulations with the 
water surface arc shown in Figure 7 and Figure 8, 
respectively. Because of the low reflectance of the 
water surface, the signal-to-noise ratios shown in 
Figure 9 are relatively low, In the visible region, the 
NESR is dominated by the photon noise, while the 
detector and electronics noise are prevalent in SWIR, 
Although we have assumed that optical efficiency 


and noise floor of the simulated sensor are 
independent of the radiation wavelength, this 
illustrates effects of spectral properties of solar 
illumination on signal-to-noise ratios in hyperspectral 
remote sensing imagery. Future simulations, 
conducted also on a subsystem level to include the 
wavelength dependence, may provide more insights 
into this problem. 
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Figure 7. Spectral radiance measured by 
the simulated satellite sensor for water . 
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Figure 8. NESR obtained for the simulated 
satellite sensor observing a water surface. 
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Figure 9, S/N obtained for the simulated 
satellite sensor observing a water surface . 
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Figure 1 L Signal-to-noise ratios obtained 
for the simulated satellite sensor observing 
dune sand. 


Mineral surfaces provide much higher signal-to-noise 
ratios than water. These are the highest S/N val ues 
among all the studied materials* S/N for calcite is 
shown as an example in Figure 10, but for all the 
studied minerals the values are similar. S/N for dune 
sand is similar to that of minerals (see Figure 1 1 ). 
Other soils display substantially lower S/N, as shown 
in Figure 1 2 for gray day as a typical example. 
Surface of concrete offers S/N, which is similar to 
that of soils, although it is slightly better in the visible 
region* Chlorophyll edge is clearly displayed in the 
S/N spectrum obtained with pine trees and shown in 
Figure 14, which is typical for almost all the 
vegetation, with exception of the dry long grass also 
studied in this group, For the man-made materials 
such as asphalt (Figure 15) and steel (Figure 16) as 
well as aluminum and terra cotta, S/N is also 
influenced by the characteristic reflectance spectra of 
the surfaces. 
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Figure 10. Signal-to-noise ratios obtained 
for the simulated satellite sensor observing 
calcite. 
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Figure 12 * Signal-to-noise ratios obtained 
for the simulated satellite sensor observing 
gray clay * 
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Figure 13. Signal-to-noise ratios obtained 
for the simulated satellite sensor observing 
concrete. 
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Figure 14, Signal-to- noise ratios obtained 
for the simulated satellite sensor observing 
pine trees. 
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Figure 15. SignaFto- noise ratios obtained 
for the simulated satellite sensor observing 
asphalt 
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Figure 16. Signal-to- noise ratios obtained 
for the simulated satellite sensor observing 
steel. 


VIL CONCLUSIONS 

We have described physics-based, end-io-end 
simulations of a remote sensing satellite system, 
which includes a hyperspectra! electro -optical sensor 
working in the visible and SWIR spectral regions. 

The simulations start with modeling of the Earth's 
surface: creation of a synthetic scene with spectral 
radiance provided for each of its elements* That 
process also includes modeling of radiance 
propagation through the atmosphere* In the next step, 
the at-sensor radiance is transformed using sensor 
effects defined by spatial sampling* optical transfer 
functions, and noise model. The calculations are 
performed for each spectral band of the sensor, so at 
the end, an image cube generated by the simulated 
sensor is created. The simulated image cubes are 
further analyzed according to a planned application, 
such as the signal-to- noise studies presented in this 
paper. These studies have shown that a hyperspectra l 
satellite sensor can provide reasonable signal -to-noise 
ratios for minerals, soils, vegetation, and man-made 
materials, while values of S/N arc crucial in accurate 
spectral un mixing and classification of remote 
sensing data. 
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